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Abstract 

We report on theoretical and experimental investigations of a novel hysteresis effect that has 
been observed on the magnetoresistance of quantum-Hall bilayer systems. Extending to these 
system a recent approach, based on the Thomas-Fermi-Poisson nonlinear screening theory and a 
local conductivity model, we are able to explain the hysteresis as being due to screening effects 
such as the formation of "incompressible strips" , which hinder the electron density in a layer within 
the quantum Hall regime to reach its equilibrium distribution. 
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I. INTRODUCTION 

decent scanning force microscope experiments and subsequent theoretical work 

3, 4 1 have shown that screening effects, notably the occurrence of "incompressible strips" 
5j, |y] , are very important for the understanding of the Hall- field and current distribution 
as well as the high precision of the low-temperature integer quantized Hall (QH) effect in 
narrow Hall samples. Since Coulomb interactions should become important within and be- 
tween the layers of an electron bilayer system showing the drag effect 

aaa 

at very low 

temperatures, we extended the theory of Ref. .4] to the bilayer case and calculated the Hall 
and longitudinal resistances for density-matched and mismatched systems. Magnetoresis- 
tance (MR) measurements on separately contacted bilayers, which were performed to check 
the results, showed pronounced hysteresis effects, similar to results reported previously for 



hole and electron double layers. In contrast to this previous work, our approach 
makes the origin of the hysteresis rather evident: the peculiar nonlinear screening properties 
leading to the occurrence of "incompressible strips" in the QH regimes of the individual 
layers. 

II. THE EXPERIMENTS 

We measured the MRs of both layers as functions of the applied perpendicular magnetic 
field and the sweep direction for matched and mismatched densities at a fixed base tem- 
perature (T = 270mK) within the linear response regime (/ ~ 50nA). We also performed 
"equilibrium" measurements, where the system was heated up to ~10K and cooled down 
again at each magnetic sweep step, in the magnetic field interval where hysteresis is ob- 
served. The samples are 15 nm wide GaAs/AlGaAs double quantum well structures grown 
by MBE and capped by top and bottom gates that control the electron densities of the 
layers. Two 2DES are confined by 5— doping remote Si donors and are separated by a 
spacer of thickness h = 12 nm. For such a separation the bilayer system is electronically 
decoupled, i.e. electron tunneling between the layers (Rtunnei > lOOMfi) is not possible and 
the system can be described by two different electrochemical potentials. Separate Ohmic 



contacts to the two layers are realized by a selective depletion technique 



12|. The samples 



were processed into 80 /zm wide and 880 /iin long Hall bars, with grown densities in the range 
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FIG. 1: Measured Hall (upper panel) and longitudinal (lower panel) resistances for the top layer 
and the bottom layer (long-dashed lines, only for the mismatched case and one sweep direction). 
No hysteresis is observed for the matched case (green dashed-dotted line) and the "equilibrium 
measurement" (solid line with circles), with n T /n B = 0.84. 

1.5 - 2.5 x 10 15 m- 2 , the mobility is 100m 2 /Vs per layer. 

Figure [T] shows Hall and longitudinal resistances of the top layer versus magnetic field 
strength, measured under different conditions. As a reference, the resistances of the bottom 
layer (long-dashed lines) are also shown for one sweep direction in the mismatched case. In 
the "equilibrium measurements" (solid lines with symbols) explained above and also in the 
case of matched densities the resistances are insensitive to the sweep direction. Note that 
the QH plateau for matched densities is narrower than that of the mismatched system in 
the "equilibrium measurement". For the density mismatched (and non-equilibrium) case, 
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FIG. 2: Longitudinal resistances versus magnetic field, for both layers and two sweep directions 
with n T Jn^ = 0.82. 

the data were taken at a sweep rate 0.01 T/min and the base temperature is always kept at 
270 mK. Apparently the resistances of the top layer follow different traces for the up- (red 
dotted lines) and the down- (black solid lines) sweep. 

Figure 121 shows that, for density- mismatched bilayers, hysteresis occurs in both layers 
under comparable conditions. The hysteresis shown in Fig. ^ occurs in the top layer near 
the high-5 edge of the v — 2 QH plateau (plateau region 5.5 T< B < 6.7 T). It occurs in 
a B-interval well within the v = 2 QH plateau (6.8 T< B < 7.9 T) of the bottom layer. 
Similarly, the MR of the bottom layer shows hysteresis at the \ow-B edge of its v — 2 QH 
plateau, in a 5-interval well within the v = 2 QH plateau of the top layer. No hysteresis 
is observed in magnetic field intervals in which the other layer is in the normal state, with 
finite longitudinal resistance. A less pronounced repetition of these features is observed in 
the v = 4 plateau regimes. 
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III. THE MODEL 



We model the electron bilayer system as a series of parallel charged planes, perpendicular 
to the ^-direction, translation invariant in the ^/-direction, and confined to \x\ < d. The 
bottom and top 2DESs with number densities n B (x) and n^(x) lie in the planes z = z B = 
and z = zt = h, respectively. Ionized donor layers with number densities and n B are 
assumed at z = — c and z = h + c, and a top gate at z = h + c + /, allowing for a density 
mismatch even with = n B = n , is simulated by ionized donors with number density 
rig(x) = n° g cosh.(^(x/d)) at z = c + h + f = z g . Solving Poisson's equation with the 
boundary condition V(x = ±d,y,z) = 0, we obtain from the charge densities ^erij(x) in 
the plane Z — Zj cLS contribution to the potential energy of an electron at position (x,y,z): 

2e 2 f d 

V 1 (x,z) = ±— / dx' K(x,x' , z, zj) rij(x'), (1) 
« J-d 

with the kernel [f| K(x, x', z, Zj) = — ln([c 2 + 7 2 ]/[s 2 + 7 2 ]), where c = cos(7r[a; + x']/4d), 
s = sin(7r[x — x']/4d), and 7 = sinh(7r[z — Zj]/4d). We write the total potential energy of an 
electron as 

V(x,z) =V B (x,z) + V T (x,z), (2) 

where V B (x, z) is the sum of the potentials created by bottom electron and donor layer, 
while V T (x, z) is the sum of the potentials due to the top electron, donor, and gate layers. 
The electron number densities in the layers are, within the Thomas-Fermi approximation 
(TEA), 

= J dED(E)f([E + V(x, Zj )-ii*}/k B T) (3) 

with j = B or T and D(E) the (collision-broadened) Landau density of states. This com- 
pletes the self-consistency scheme j^, Q|- 111 the practical calculations we first decoupled 
the layers, replacing V(x,zj) by V 3 (x,Zj), and solved the single electron-layer problem for 
bottom and top system separately. There it is convenient to fix the edges x = ±6 of the 
electron density profile in the limit of zero magnetic field and temperature j^lQ], which fixes 
the average density (and n BT ). With the converged results at finite T and B, we treat the 
inter-layer coupling iteratively. 
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To describe the density-mismatched case, we add more electrons to the top layer by 
setting Vo = n®/n to a positive value and keeping the depletion length d — b fixed. We scale 
energies by the average Fermi energy E* F = (Ep + E§)/2 , e.g. tkjj c /E* F = Q C /E F . 

Next we implement in each electron layer a quasi- local model to impose a current, fol- 
lowing the lines of Ref. |4j. It is based on a local version of Ohm's law, j(x) = a(x)E(x) 
between current density j(x) to driving electric field E(x). The position-dependent con- 
ductivity tensor a(x) is taken from a bulk calculation 13], with the electron density n e \ 
replaced by the local one, and a spatial smoothening over a length of the order of the Fermi 
wave length is performed, in order to avoid unphysical artifacts of a strictly local transport 
model and of the TFA [4]. Typical results with a slightly modulated donor distribution 
tIq{x) = tIq(x) = n [l — 0.01 cos(57nr/2cZ)], introduced to simulate the long-range-part of 
the impurity potentials |4j , are shown in Fig. El 

IV. SIMULATION OF NON-EQUILIBRIUM 

The quasi-equilibrium model just introduced is, of course, not able to describe hysteresis 
effects. It yields, however, qualitatively correct results for the matched density case (n° = 0), 
where the MRs for both layers are identical, and for the mismatched systems it yields results 
in qualitative agreement with those found in the "equilibrium measurements" . Moreover, the 
results shown in Fig. El give some hints towards the possible origin of the hysteresis effects. 
We see that, in the plateau regime of the QH effect, i.e., in the magnetic field interval in 



which incompressible strips (ISs) exist J], the position of the incompressible strips [local 
filling factor u(x) = 2] and the potential distribution change drastically. To realize these 
changes, electrons must be transported, e.g. by thermal activation, across the ISs. In real 
samples (of finite length) this is extremely difficult, since states at the Fermi level exist on 
both sides of an IS but not within the IS, and since the region between ISs is not connected to 
the contacts. Therefore, after a slight change of the magnetic field, it may take an extremely 
long time until the electron density relaxes to its new equilibrium distribution . 

To simulate this hindered approach to equilibrium, we calculate the MRs of the current- 
carrying layer in a 5-interval in which the other layer exhibits a QH plateau, using different 
frozen-in density and potential distributions for up- and down-sweeps of the magnetic field 
B. For the down-sweep, we freeze in the density profile of the other layer at the high-5 
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FIG. 3: Density profiles (color scale) across the sample as a function of magnetic field (left panel) 
and (normalized) Hall potential profile of both layers for selected B values (right panel). The 
density mismatch is obtained by setting Vq/Eq = 0.05 which results in nfjn^ = 0.84 at a low 
temperature ksT/Ep = 0.0124 for fixed b/d = 0.9. Here Eq(= 2Tre 2 nod/R) is the pinch-off energy 
which defines the minimum of the bare confinement potential. 

edge of its QH plateau and for the up-sweep we freeze it in at the \ow-B edge, and use these 
fixed potentials to describe the other layer when calculating the equilibrium and transport 
properties of the current-carrying layer. Thus, for the case shown in Fig. |2l we calculate the 
MRs of the bottom layer in the QH plateau regime of the top layer by fixing the potential 
due to the top layer at Q c / ' E* F w 1.3 for the down-sweep and at Q C /E F w 1.0 for the up- 
sweep. To calculate the MRs of the top layer in the plateau regime of the bottom layer, 
we fix the potential profile of the bottom layer at Q C /E F ps 1.16 for the down-sweep and 
at Vt c /E* F ps 1.05 for the up-sweep. Since the MRs slightly outside the QH plateau regimes 
depend on the electron and current density profiles [4], we will obtain slightly different 
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FIG. 4: Longitudinal resistances for high density mismatch parameter. The hysteresis-like effect 
is observed in the active layer only if the passive layer is within the plateau regime. 

results for the different sweep directions. Figure 0] shows longitudinal magnetoresistances 
calculated along such lines. 

V. SUMMARY 

We report on measurements and calculations of the magnetoresistances in separately 
contacted e-e bilayer systems. The calculation, based on a self-consistent Thomas- Fermi- 
Poisson theory of the equilibrium state and a quasi- local transport theory Q, yields rea- 
sonable agreement with the measurements for density-matched layers and for mismatched 
layers, if the measurement allows for thermal relaxation at each value of the magnetic field 
B. Z?-sweeps at constantly low temperature show hysteresis between up and down sweep 
in the current-carrying layer, provided the other layer exhibits the quantized Hall effect. A 
model, based on the existence of incompressible regions in the quantum Hall states j^J and 
the extremely long relaxation times in compressible regions surrounded by incompressible 
ones is worked out and can explain the observed hysteresis effects. 
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